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J>^! Abstract 
■ 

We consider brane world models with interbrane separation stabilized by the Goldberger- 
Wise scalar field. For arbitrary background, or vacuum configurations of the gravitational 
and scalar fields in such models, we construct the second variation Lagrangian, study its 
gauge invariance, find the corresponding equations of motion and decouple them in a 
. suitable gauge. We also derive an effective four-dimensional Lagrangian for such models, 

1^ I which describes the massless graviton, a tower of massive gravitons and a tower of massive 

' scalars. It is shown that for a special choice of the background solution the masses of the 

graviton excitations may be of the order of a few TeV, the radion mass of the order of 100 
GeV, the inverse size of the extra dimension being tens of GeV. In this case the coupling 
l/^ I of the radion to matter on the negative tension brane is approximately the same as in the 

■ unstabilized model with the same values of the fundamental five-dimensional energy scale 

. and the interbrane distance. 

^ ; 1 Introduction 
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Brane world models and their phenomenology have been widely discussed in the last years 
[T]-[6]. One of the most interesting brane world models is the Randall- Sundrum model with 
5^ i two branes, - the RSI model [7]. This model solves the hierarchy problem due to the warp 
factor in the metric and predicts an interesting new physics in the TeV range of energies. A 
flaw of the RSI model is the presence of a massless scalar mode, called the radion, which arises 
due to the fluctuations of the branes with respect to each other. Its interactions contradict 
the existing experimental data, and in order the model be phenomenologically acceptable the 
radion must acquire a mass, which is equivalent to the stabilization of the brane separation 
distance. The latter can be achieved by introducing a five-dimensional scalar field with bulk 
and brane potentials, whose vacuum energy has a minimum for a certain interbrane distance 
[8] . A disadvantage of the approach proposed in [8] is that the backreaction of the scalar field on 
the background metric is not taken into account. This problem is solved in the model proposed 
in [9]. 

Nevertheless, most of the papers on the phenomenology of the RSI model consider the 
unstabilized model, just putting the radion mass by hand. Such an approach seems to be 
inconsistent, because the backreaction of the scalar field on the metric leads to a renormalization 
of the parameters of the RSI model. The scalar sector of the stabilized RSI model in the 
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background of [9] was studied in [TO], [TT], where the fundamental five-dimensional energy scale 
of the theory was assumed to be of the order of four- dimensional Planck mass. In [12] it was 
shown that a consistent physical interpretation of the theory on the negative tension brane for 
the model of [9] is possible, if the fundamental five-dimensional scale is of the order of ITeV, 
rather than the Planck one. 

In the present paper we are going to study all the physical degrees of freedom of the stabilized 
RSI model, i.e. both the tensor and the scalar sectors, in an arbitrary background. Our 
approach is based on the Lagrangian description of the linearized gravity, which was developed 
for the unstabihzed model in [13]. We find a convenient gauge and decouple the equations for 
tensor and scalar modes. Then for the background solution of [9] we calculate the masses of 
tensor and scalar excitations and their couplings to matter on the negative tension brane. 



2 Linearized gravity in stabilized brane world models 

Let us denote the coordinates in five- dimensional space-time E = x S^/Z2 by {x*^} = 
{x^, ?/}, M = 0, 1, 2, 3, 4, /i = 0, 1, 2, 3, the coordinate x^ = y, —L < y < L parametrizing the 
fifth dimension. It forms the orbifold, which is realized as the circle of circumference 2L with 
the points y and —y identified. Correspondingly, the metric qmn and the scalar field satisfy 
the orbifold symmetry conditions 

9t^uix, -y) = g^,u{x, y), g^^ix, -y) = -g^A^x, y), 

guix, -y) = gu{x, y), -y) = y). (1) 

The branes are located at the fixed points of the orbifold, y = and y = L. 
The action of stabilized brane world models can be written as 



where Sg and Ss are given by 



S — Sg + S^, (2) 



Sg = 2M-' I d'^x I dyRy/^, 



S^ = - I d^x I ^dy Q(/*^^aM05^0 + ¥{(!)) ) 



g\i{(j))d'^x — I ^/—gX2{4')d'^x. (3) 
y=o Jy=L 

Here V{(j)) is a bulk scalar field potential and Xi^2{4>) are brane scalar field potentials, g = detg^^, 
and cjfj^u denotes the metric induced on the branes. The signature of the metric guN is chosen 
to be (-,+,+,+,+). 

The standard ansatz for the metric and the scalar field, which preserves the Poincare in- 
variance in any four- dimensional subspace y = const, looks like 

ds^ = e-^^^y^r^^^dx^'dx^ + dy"^ = -fMN{y)dx^^dx^ , (4) 
(f){x,y) = (j){y), (5) 
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T]f^i, denoting the flat Minkowski metric. If one substitutes this ansatz into the equations 
corresponding to action ([2]), one gets a rather comphcated system of nonhnear differential 
equations for functions A{y), (f){y): 

i2M%Ar+^{v-^i<pr) = o, 

Q(0O' + V + Xi6{y) + X26iy - L)^ = -2M'' {-3A" + 6{A'f) . (6) 
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Here ' = 84 = d/dy. 

Suppose we have a solution A{y),(l){y) to this system for an appropriate choice of the 
parameters of the potentials such that the interbrane distance is stabilized and is equal to L. It 
means that the vacuum energy of the scalar fleld has a minimum for this value of the interbrane 
distance. 

Now the linearized theory is obtained by representing the metric and the scalar field as 

9MN{x,y) = ^MN{y) + ^y==hMN{x,y), (7) 

^{x,y) = (f){y) + -jl=f{x,y), (8) 

substituting this representation into action ([2]) and keeping the terms of the second order in 
Hmn and /. The Lagrangian of this action is called the second variation Lagrangian and has 
the form 



-2 Va/ h''^ hsN - Vsh h) + {A')' HhMNh^''' - hh) - 
-A" (hMNh'''' - fhh + '^hMuh'^') + 4173 {hMNh^''' - \hh) + 
+ \ [Kuh^^- - \hh) [\^6{y) + \26{y - + 

-d'^fduf -P{^ + ^Siy) + p6iy - L))]. (9) 

Here h = 'jMNh^^'^ , h = •y^uh'^", (p stands for the background solution and Vat denotes the 
covariant derivative with respect to metric •Jmn- 

Varying the action built with this Lagrangian and taking into account background field 
equations ([6]), we arrive at the following equations of motion for the fiuctuations of metric: 

1. /iz/-component 
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+\l^,v(y-^A^Jl - d^d'^hu - d^d'^h + AAdJi - 3A'94/i44+ 
+|/i447m.(4^" - - (/>' + /(-4AV' + 0")) = 0; (10) 

2. /i4-component 

+2]^^^/'/'' = 0; (11) 

3. 44- component 

d^'id'^h^^ - (9^/i) - QA'd>'h^i + 3A'd^h+ 

+ 2W [-h^^y -f% + /V') = 0; (12) 

4. equation for the field / 

hAA [^5{y) + ^5(y - L)) - 2/^44(-4A>' + 0")+ 
+(/.'a4^ - 0'a4/i44 - 2(j)'d''h^4 - SA'dJ + 2dMd^^f- 

-2/ (0 + + pSiy - L)) = 0. (13) 

We will also use the following auxiliary equation, which is obtained by contracting the indices 
in the /ii/-equation: 

d^d^hf,^ - d^'d^ih + - QA'd^^Ki -h)- Id^dji - 6A'(9^/i^4+ 
+?>d^^d,h^, + 6/i44(4A'2 - A") - ^(/V + /(-4A'0' + 0")) = 0. (14) 

These equations were also discussed in [12]. The normalization of fields in the above equations, 
as well as in Eqs. ([7]), ([H]) and (jH]) differs from the one adopted in this paper. It is more 
convenient in the Lagrangian approach and does not effect the equations of motion. 
These equations are invariant under the gauge transformations 

hll{x, y) = hf,^{x, y) - {df,C.u + d^i^, - 27^^(94A^4), 
y) = h^i{x, y) - {d^^^ + d^^^, + 2diA^^), 

h^ilix, y) = /i44(x, y) - 284^4, (15) 
/W(x,y) = /(x,y)-940e4, (16) 

provided C,M{x,y) satisfy the orbifold symmetry conditions 

^,i{x,-y) =^f,{x,y), U{x,-y) = -U{x,y). 
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These gauge transformations are a generalization of the gauge transformations in the unsta- 
bihzed RSI model [61 [ISj. We will use them to isolate the physical degrees of freedom of the 
fields Hmn and /. Let us show that the gauge transformations with function ^4 allow one to 
impose the gauge condition 

{e-'^h,,y - ^,e-'Wf = 0. (17) 
Really, Eqs. ffT5]) . ffTBl) imply the following equation for ,^4 

94^4^4 - 2A'd,^, - = ~{d,hu - 2A'hu - ^07). 

The standard theory of differential equations demands the coefficient in front of ^4^4 to be 
continuous, though in our case it is not so. To cure this flaw, we rewrite the equation as 
follows: 

d.d.i^e'^) + i,e~^ (a" - {Af ^ <-'^^ 



6M3 

= -^6-^(^4/144 - 2A'/i44 - ^</>7) = w{x,y). (18) 

The functions {A'Y and (0')^ are smooth functions of y. Although function A" has 5-like 
singularities at y = 0,y = L, the singular terms drop from the equation, because ^4 is equal to 
zero at these points. Thus, the factor in front of ^4e~^ is a smooth function, and the equation 
for 1^46-"^ in the interval [—L, L] can be treated by the standard methods. In accordance with 
the general theory, the homogeneous equation, corresponding to f|T8|) . has two independent 
solutions in the interval [0,L]; we denote them Xiiv) a^id X2{y)- Now we can use the Green 
function method to find ^4, which gives: 

Jo W{xi{z)X2{z)) 

A xi{y)x2{0) - X2{y)xii0) xi{l)x2{z) - x2{l)xi{z) 



Xi{L)x2{0) - X2{L)xm Jo W{xi{z)x2{z)) ^''^ 

Here W{xi{z)x2{z)) stands for the Wronskian of the solutions Xiiv) and X2{y), which in the 
case of Eq. f|T8l) is just a constant. It is not difficult to check that the obtained function ^4 is 
equal to zero at the ends of the interval and therefore can be continued to an odd function on 
[— L, L]. Thus, gauge condition ( ITTI) really exists. We also note that this relation was obtained 
in [To] from the equation for /i4-component, in which only the scalar degrees of freedom were 
retained. Similar to the case of the unstabilized RSI model, the gauge transformations with 
functions allow one to impose the gauge h^i{x,y) = 0, after which there remain the gauge 
transformations satisfying 

d,ie'%) = 0. (20) 

Thus, we can use the gauge 

(e-^^/^44)' - = 0, 

V = 0. (21) 



5 



Next we represent the gravitational field as 

hf,u = V + ^7m^^^, (22) 

with bfj,i, being a traceless tensor field {^'^"h^y = 0). 

Substituting gauge conditions and representation (1221) into the /i4-equation and into 
contracted /iz/-equation f[T^ . we get: 

- 94(5^6^.) + ^dMh + 2/^44) = 0, (23) 



d'-d'^b,^ - -d,d^h - -d,d^h,, - -—h + 

+QA'dJi - 3—^44 + I'lA'dihii = 0. (24) 



Eq. ( l23l) suggest the substitution /i = —2/144, which allows one to decouple the equations 
for the fields b^^, and /. The possibility of using this substitution and representation (!22|) 
for decoupling tensor and scalar equations was mentioned in [llj . 

Really, as a result of this substitution Eqs. (!23l) . take the form 

d.id-'b^.) = 0, (25) 
d'^d'^b^, = 0. (26) 

It is not difficult to check that the residual gauge transformations fl20l) are sufficient to 
impose the gauge |13] 

d'b,, = 0, 

in which the former equations are satisfied identically. 
Thus, in what follows we will be working in the gauge 

(e-^^/^44)' - = 0, 

= 0, 

b = r^b^. = 0, d^b^, = 0, (27) 

the residual gauge transformations now being 

= e-^\ix), d'e.ix) = 0, De, = 0. (28) 

Obviously, after the substitution h = —2/144 contracted /iz/-equation (IT^ and the /i4- 
equation are satisfied identically in this gauge. Eq. (ITUI) for the /iz/-component reduces to 
an equation for a transverse-traceless tensor field b^u{x,y): 



i(e-(^)nV + ^)-&..(2(A')^-^")=0. 



(29) 



This equation does not include the background scalar filed (j){y) and is absolutely analogous to 
the corresponding equation in the unstabilized RSI model. 
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In order to find equations for the scalar field h^^, we have to solve gauge condition (l27Il with 
respect to / and to substitute the latter into Eq. (IT^ . (IT^ . This can be done either using the 
regularization {sign{y)Y = 1, or restricting the equations to the interval (0, L) and taking into 
account their singular terms with the help of boundary conditions. The latter technique turns 
out to be simpler and we will use it. 

Equation for 44-component (|T2l) simplifies considerably, when rewritten in the interval (0, L) 
in terms of a new function g = e~'^'^^^^h44(x,y) and using the expression for the potential V in 
terms of A and d) (EI): 



(30) 



Let us note that substitution (!22|) and gauge condition (1271) . which decouple the equations of 
motion, in terms of g take the form: 

h^,u = V - ^V^Lug, hu = e^^^y^g, (31) 
' ^ e-2V/ = 0, (32) 



" 3M3 

^M = 0, 6 = 7m-&''' = 0, d''h^, = Q. (33) 



Gauge condition fl32l) solved for / in the interval (0,L), looks like 



/ = 3M^^^'. 

Substituting this expression for / into Eq. f|T3l) gives an equation, which is obtained by differ- 
entiating Eq. ( !30l) with respect to and boundary conditions on the branes: 



2A 



g' + d,d''g\y=+o = 0, 



Id'X 



2 



/>2 



+ ^]g' -d,d^'g\y=L-o = 0. (34) 



Thus, our gauge choice and the substitution enabled us to decouple the equations of motion. 
This means that the fluctuations of the metric and of the scalar field against any background 
given by a solution to Eqs. ([6]) are described by two fields, - tensor field b^^{x, y) and scalar field 
g{x,y). Their classical equations of motion are given by Eqs. (12^ and (15(1 . respectively. 

3 Mode decompositions 

Let us study first the modes of the tensor field 6^i,(x, y), which satisfies Eq. ( !29|) . Substituting 
into this equation 

b^uix,y) = c^^e'P'^iJniy), c^^ = const, p"^ = -ml, 
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restricting it to the interval (0, L) and replacing the singular terms by the boundary conditions, 
we get: 



< + 2A'tPn\y=+0 = < + 2A>„|^=L_o = 0. (35) 

The boundary conditions suggest a substitution i/jn = exp{—2A)ujn, which turns this equation 
into 



^'n\y=+0 = i^n\y=L-0 = 0. (36) 



We see that the eigenfunctions ujn are solutions of a Sturm-Liouville problem with von Neumann 
boundary conditions. In accordance with the general theory, |T3] the problem at hand has no 
negative eigenvalues for arbitrary A, only one zero eigenvalue, corresponding to ujq = const, 
and an infinite number of positive eigenvalues, asymptotically given by the formula 

ml = ^, l = j\^(yUy. (37) 

This formula should be specified for finding the masses on different branes. We recall that the 
masses of excitations on each brane should be calculated in the Galilean coordinates [6l [121 [13] ) 
for which A{y) is equal to zero on the corresponding brane (we recall that coordinates are 
called Galilean if g^^, = (—1, 1, 1, 1) [15]). Thus, the latter formula can be explicitly adopted 
for calculating masses on the branes as follows: 

/ = e^^^y^-^'^^^^dy for the brane at y = 0, (38) 
/ = /^^e(^(^)-^(^))d|/ for the brane at y = L, (39) 

which is valid for an arbitrary A{y) satisfying Eqs. ([6]), because they define it up to an additive 
constant. 

The eigenfunctions {i'n{y)} of eigenvalue problem (1531) build a complete orthonormal set, 
the eigenfunction of the zero mode being 

My) = Ne-'^^y\ (40) 

Expanding b^^ in this system 

oo 

K^ = Y.K'^^^^My), (41) 

n=0 

we get four-dimensional tensor fields b''^^{x) with definite masses. An important point is that 
due to the form of the zero mode eigenfunction residual gauge transformations f[28|) act only on 
the massless field b^^j^{x) and provide the correct number of degrees of freedom of the massless 
graviton 
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In order to find the mass spectrum of tfie scalar particles described by Eq. (130|) let us 
substitute 

g{x,y) = e'P''gn{y), = -fij 



into this equation. As a result, the equation and the boundary conditions for gn{y) take the 
form: 

g: + 2Ag: - 2'^g: - = -/x^e^^., (42) 



6M3 

1 <P\i 
2 



2 -lj]9'n + fiy^gn\y=+0 = 0, (43) 



' +^]g'r,-fiy^gn\y=L~0 = 0. (44) 



Let us write Eq. (I42l) in the Sturm-Liouville form: 



;9n] - 7^9n = (45) 



It is not difficult to see that the operator in the eigenvalue problem for this equation with 
boundary conditions ( H3i) . ( Hil) is not self-adjoint. Nevertheless, these boundary conditions look 
like the usual Sturm-Liouville boundary conditions and lead to a number of general assertions 
about the spectrum and the eigenf unctions of the problem at hand. 

Multiplying Eq. (H5|) by (jn, then integrating over (0, L) and integrating by parts in the 
term with derivatives, we get: 



. , , ,2^ A^^'Ai 0"V' , ,2, 

1 /_ 7Ji\^^9n\ dy+ (^2^ JIi\^\9n\ \y=+0+ 

ld^\2 , (t)"Y' e^^ 
2 



+ 1 o J^2 + ^/ /A'\2\9n\ \y=L-0 



1 



6M3 Jq 



e^^\gn?dy+ / -7--;r^\g'n?dy. 



This means that if 



all the eigenvalues of the eigenproblem are real and positive. The standard technique for proving 
the orthogonality of eigenf unctions gives in this case for m ^ n 

io (0')2^'"^-^^ - ^2 <P' ) ^^'y9mgn\y=+0 + 
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Thus, the eigenfunctions of this problem, corresponding to different eigenvalues, are not or- 
thogonal with respect to the weight suggested by the form of Eq. (HSj) . 

It is very difficult to prove rigorously that this set of eigenfunctions is complete. We can just 
argue that due to the special form of boundary conditions (H3|) . fHil) the set of eigenfunctions of 
the eigenvalue problem under consideration is complete. Really, for the parameters of the scalar 
field potential ^—^^ oo the eigenvalues drop from the boundary conditions, and the operator 
in the equation becomes self-adjoint (this approximation was used in [TU], [H]). Therefore, for 
— > oo the eigenfunctions of the problem under consideration go to the eigenfunctions of 
a Sturm-Liouville problem with a self-adjoint operator, which build a complete orthogonal set. 
The orthogonality of the eigenfunctions for | '^^^2 ^ 00 can be also seen in Eq. (H71) . Thus, we 
assume that the eigenfunctions {gn{y)} of the problem (142|) - fHi|) form a complete denumerable 
non-orthogonal set. We would also like to note that for /i,^ 00 boundary conditions (HSl) . 
(jH]) go to 5'„(0) = gn{L) = 0, and therefore for large m or n the integral in (H7I) is close to zero. 
Such a set of functions may be called asymptotically orthogonal. 

It is easy to understand that the eigenfunctions Qniv) can be chosen to be real. Then 
five-dimensional scalar field g{x, y) can be expanded in this system as 



n=l 



where the four dimensional real scalar fields (pn{x) have masses /x^. 

Now we can find the effective four-dimensional action for the system. Substituting fl22l) . 
( l3Tl) -( !33|) into Lagrangian (Q, we find that the tensor field and the scalar field Lagrangians 
decouple. Then substituting expansions f HTl) . f HHl) into these Lagrangians and using Eqs. ( l35l) . 
( H2l) . assuming the eigenfunctions {ipn} and {gn} to be appropriately normalized and integrating 
over the coordinate of the extra dimension, we find that the reduced tensor field Lagrangian is 
a sum of the standard four-dimensional Lagrangians of tensor fields with masses m„, whereas 
for the scalar field Lagrangian we get 

^scalar = [v'"' ^l^^findu^k + f4^n^k] dyC^^ (^^gk + J^gWk 

Now consider the integral over dy. Integrating by parts the term with derivatives, using Eq. 
[5]) and boundary conditions (H3l) . (jH]), we arrive at the result 



7Ti^^9ngk\y=+o+ 



A2 . e^^ 



^ ' 2 ^ ' ( M\23n9k\y=L~0 

For real gn{y) Eq. (H7|) now implies that this expression is equal to zero for n ^ k. Assuming 
the eigenfunctions {gn\ to be appropriately normalized, we find that the reduced scalar field 
Lagrangian is also the standard Lagrangian, and the complete reduced action is 
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— J dx {du^kd''ipk + fJ'l^k'fk) ■ (49) 

k 

We see that the effective Lagrangian of the tensor fields coincides with the one of the unstabihzed 
model, [13] whereas the effective Lagrangian of the scalar fields coincides with the one found in 



TT] for the case ^-n^ oo. 



Now we are able to find the couplings of the four-dimensional fields b'^^^^x) and fn{x), 
to matter on the branes, which are defined by the coupling of the fluctuations of the five- 
dimensional gravitational field /i^i/ to matter on the branes. The latter is given by 



L= [ h^^{x,Q)Tl^J^-det^,Mdx + 



1 

+ , 



h^^^ix, L)T^f J-det^^,{L)dx, (50) 

B2 ^ 



where T^^ and are energy-momentum tensors on brane 1 and 2, respectively and 



y(l,2) 



Qg^U ^^^^ 



for our choice of the metric signature. 

In what follows we restrict ourselves to considering the fields on the brane at y = L only, 
which we assume to be "our" brane. Substituting decompositions fHTl) . fllSl) into fISUl) we find 
that the interaction of the tensor and the scalar fields with matter on the brane at y = L in 
the Galilean coordinates is 



V8M3 



B2 



n=l 

oo \ 

--Y,9n{L)Mx)Ti:]dx. (51) 



2 

n=l 



Thus, the couplings are defined by the values of the wave functions on the brane. The latter 
can be found only if we specify the background solutions A{y) and (p{y). 



4 Specific example 

The considerations in the previous sections allowed us to find the general structure of the brane 
world models stabilized by the scalar field. To make predictions for the masses and the coupling 
constants we must take specific potentials for the scalar field and a particular vacuum solution 
of system IQ. 

To find an analytic solution to this system we will use the results of [9], [16]. Let us consider 
a special class of potentials, which can be represented as 
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It is easy to check that if we put 

1 dW 1 
(p'{y) = -^^aniy)-^, A'{y) = sign{y)^^W{(l)), (52) 

then Eqs. ([6]) are vahd everywhere, except for the branes. In order the equations of motion be 
valid everywhere, one needs to finetune the brane potentials Xi,2{(p)- 
Let us take W^cj)) to be 

W = 24:M^k - (53) 
so that V{(j)) is a quartic potential. Finetuned potentials on the branes can be chosen as follows: 

Ai = W^(0i) + W'{<PM - 0i) + /32(0 - 0i)2, (54) 

A2 = -1^(02) - H^'(02)(0 - 02) + f3ii(l) - 02)'. (55) 

The parameters of the potentials /c, m, 0i,2, /5i,2, when made dimensionless by the fundamental 
five-dimensional energy scale of the theory M, should be positive quantities of the order 0(1), 
i.e. there should be no hierarchical difference in the parameters. 

For such a choice of the potentials the solution of the equations of motion is given by [9J 

0(y) = 0ie-"IJ^I (56) 

The interbrane distance is defined by the boundary conditions for the field and is expressed 
in terms of the parameters of the model by the relation 



llnf^V (58) 



U \(P2 

Thus, we see that the brane separation distance is stabilized. 

Let us study the mass spectrum of tensor particles, which is defined by Eq. fl35l) . The zero 
mode solution for arbitrary A is given by fBUl) . For our choice of A it is impossible to find exact 
solutions for other modes. Therefore, in what follows we will use an approximation uL <^ 1, 
which is rather general and physically interesting [12j. Keeping in A only the terms linear in 
y, we get ^ 

A{y) = ~k\yl ~k = k-^u. (59) 

Thus, in this approximation Eq. (1291) for the tensor field coincides with the equation of the 
unstabilized model, where a substitution k ^ k was made. This equation can be solved 
exactly, and the formulas for eigenfunctions and eigenvalues were discussed in detail in [13]. In 
particular, in this approximation the normalized functions ^|Jo for mo = look like 

My) = Noe-''^y\ No= ^\ (60) 
in the coordinates, which are Galilean on the brane at y = 0, and 

My) = Noe-'^\y\^'^\ No = ^J" . (61) 
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in the coordinates, which are Gahlean on the brane at y = L. For the first massive tensor 

excitation on the brane at y = L we get [13] ipi{L) ^ —\fk, m\ ^ 3.83A;. Eqs. ( [37l) . ( !38l) . ( l39l) 
give the following mass spectra for large n: 

ml = TT^Pn^e"^^-^ (the brane at y = 0), 

ml = -k'^k^v? (the brane at y = L). (62) 

Let us study now the mass spectra of the scalar particles. Equation and boundary conditions 
( H2l) - ( l44l) for our choice of the potentials and background solutions take the form: 

9n + ^^9n- ^^9n " -^]^9n = -f^n^ 9n, (63) 

iPl + u)g'^ + fxy^g„]y=+o = 0, (64) 
(/?2 - u)g'n - file'^^gn\y=L-o = 0. (65) 

To find explicitly the eigenf unctions and eigenvalues of the problem we will use the same 
approximation uL -C 1, which now implies 

y 

uy = 'uL— < uL <ti 1. 
Substituting the explicit form of (p into fl63|) we get 

gl + 2Ag'^ + 2< - -^,u^9n + y^y9n = 0. (66) 

We solve this equation in the coordinates, which are Galilean on the brane at y = L. In 
this case A{y) is expressed in terms of k ( l59l) as: 

Aiy) = k{y-L). 

We have already found that for /9| — u > all eigenvalues of this problem are larger than zero. 
Therefore, we introduce a new variable by the relation 

k k k 

In terms of this variable the equation takes the form: 

Let us look for gn in the form gn{z) = z°'tn{z). Then the equation for t„ is 



+ ( 2a + 3 + — ) + ^ ^5 '-tn+ [l-—]tn = Q. {6 



dz^ V k J z dz z^ \ z 
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In order to turn this equation into the Bessel equation we put a = — (^1 + and get 

The general solution to this equation is 

tn{z) = AJaiz) + BJ^aiz), a = Va?~+¥ ■ 

Correspondingly, we get 

9n{z) = z-('^f) (AUz) + BJ^^iz)) . (70) 
The boundary conditions in terms of z look like: 

kz'^gn + (/?! + u)z^\^^^i^^.kL = 0, 

kz'^gn - ip2 - u)z^\^^^i^ = 0. (71) 

Below we show that for reproducing the Newtonian gravity on the brane ai y = L for strong 
five-dimensional gravity we must take kL ~ 35. In this case zi = ti^^-^^ ^ is a very good 
approximation, and the boundary condition at zero allows us to drop the singular term with 
J~a{,z) in gn{z) since B /A ~ e~^^^, and the corrections due to this term are negligible. Thus, 
up to normalization gn{z) can be written as 



9n[.Z) 



^^"^hjz). 



The second boundary condition at 2:2 = ^ gives an equation for the mass spectrum of the 
scalar particles: 



u kz. 



l + a+y + J^{Z2) - Z2Ja-l{z2) = 0. (72) 

k f3i-uj 

Expanding the Bessel function for small Z2 up to the second term and keeping the terms up 
to the order Z2 in the equation, we get the following relation for the mass of the lowest scalar 
excitation: 

2 ^ 4fc^(-l+a-f)(l + «)(/?! -tx) 
4fc(l + a) + (l + a- !)(/?!-«)■ 

For small u/kit reduces to 

u? = ^ (74) 

which for /5| ^ cxd formally coincides with the results of [TU], [TT], [12] in our approximation 
(we recall that the fundamental energy scale in [TUj, [H] is the Planck one). The next roots of 
Eq. (172!) are of the order /c^ , and the asymptotic formula of the general theory [H] for large n 
gives /i^ = n'^kP'n^. 
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The normalization condition for the eigenf unctions {gn{y),n = 0, 1, ...} of the problem (jl2 
which gives the canonical kinetic terms in fl49l) . is the following: 



dye'^'^ ( gngk + 



The normalized functions gn{y) look like 



9n9k 



(75) 



(76) 



Ar. 



U(f)i [fln/k 



_\ 1+u/k 



-uL 



X 



X 



1 

2k 



k + U ^ fin 



k^a^ 



+ 1- 



PI 



u 



We also note that in our approximation the boundary term at j/ = drops from the normal- 
ization condition. 

Now that we have found explicit expressions for the wave functions of the tensor and the 
scalar fields, we can find their coupling constants to matter on the brane at y = L. 

The coefficient in front of the zero mode of the tensor field &°i,(x) in ( l50i) can be expressed 



(77) 



in terms of the Planck mass on the brane at y = L as 1/ a/ 8Mpi . Then Eq. f l6T|) implies 



Nn 



k2 



V8M3 



-1) 



1 5 
2 



which gives a relation between the Planck mass on the negative tension brane and the funda- 
mental five- dimensional energy scale M 



Mil 



- 1 



k 



(7J 



Again we have obtained the formulas, which coincide with the formulas of the unstabilized 
model after replacing k k. Therefore, in order to get the correct value of the gravitational 
constant on the brane aX y = L one can take the same values of the parameters, which are used 
in the unstabilized model, i.e. M ~ /c ~ ITeV, kL ~ 35 [12]. Then the parameter u may be 
of the order of ten GeV, and the radion mass fl74l) may be of the order of hundred GeV. 

The exact solution for the zero mode of the tensor field f HOl) allows one to get the relation 
between the Planck mass on the brane aX y = L and the fundamental five-dimensional energy 
scale beyond the approximation used here. Other possible approximations were studied in [12] . 

Eq. (1511) shows that the coupling constant of the n-th scalar mode to matter is defined by 
the value of its wave function on the brane at y = L and is approximately given by the relation 



gn{L) 



Ar 



2V8M3 



,j I ^ 

V32AP "Kk 



/in 
k 



-{1+u/k) 
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For the radion, this constant is 



k 



ei 



24M3' 



and turns out to be of the order ^ ~ 5TeV for the above given values of the model parameters. 
One can see that ei does not depend on parameter (32- 



5 Discussion 



In the present paper we have considered the general structure of the brane world models sta- 
bilized by the scalar field. For an arbitrary background configuration of the gravitational and 
scalar fields, satisfying the equations of motion for a stabilized model, we constructed the sec- 
ond variation Lagrangian and derived the equations for the fields describing the fluctuations 
against the background. A convenient gauge and a substitution were found, which enabled us 
to decouple the equations of motions and to isolate the five-dimensional degrees of freedom. 
It was shown that the tensor sector splits from the scalar one and has the same structure, as 
in the unstabilized model. Namely, for any background there is a massless four-dimensional 
graviton and a tower of massive tensor fields, which represent the four- dimensional degrees of 
freedom in the tensor sector. 

The structure of the scalar sector was found to be more complicated: the operator of the 
mass squared turned out to be non-self-adjoint and its eigenfunctions to be non-orthogonal. 
Nevertheless, the structure of the five-dimensional action is such that the non-diagonal inter- 
action terms of the four-dimensional scalar modes vanish. 

For a particular choice of the background solution and the parameters of the model it was 
found that the influence of the scalar field background on the tensor excitations reduces to a 
renormalization of the parameter k of the unstabilized model, which is replaced by k (l59l) . In 
this case the inverse size of extra dimension may be of the order of tens of GeV and the masses 
of tensor excitations may be of the order of a few TeV, the radion mass being of the order of 
lOOGeV^. For this choice of the parameters the radion coupling constant to matter on the brane 
aX y = L turned out to be of the same order as in the unstabilized model with the same choice 
of the parameters. 

At the same time it is quite possible that it can be much larger for a different choice of the 
model parameters. In this case the radion mass and the masses of the tensor excitations must 
also shift. To find out, whether it is really so, one has to scan the whole parameter space of the 
model. The coupling of the radion to matter can also become stronger due to the radion-Higgs 
mixing as discussed, for example, in [TD], [T7]-[19j. 
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